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REFRACTED-RAY  SCANNING  (REFRACTED  NEAR-FIELD  SCANNING) 
FOR  MEASURING  INDEX  PROFILES  OF  OPTICAL  FIBERS 

M.  Young* 

National    Bureau  of  Standards 
Boulder,   Colorado  80303 

The  purpose  of  this  work   is  twofold.     First,    it  provides  an  elementary 
description  and  tutorial    overview  of   the   new  ref racted-ray  method  of  measuring 
fiber-index  profiles.      Second,    it  presents   new  results  concerning  the  theore- 
tical   foundation,    the   linearity   and  precision,    resolution   limit  and  edge 
response,   and  other  aspects  of  the  method.      In   particular,   we  find  that  index 
differences  may  be  measured  to  5  percent  or  better  and   spatial    resolution   is 
diffraction   limited.     We   conclude  by   showing  about  3  percent  agreement  with 
another  laboratory  and  good  agreement  with   numerical -aperture  measurements  per- 
formed by   participants   in  a  round-robin  experiment. 

Key  Words:  Fiber  index  profile;  index  profile;  near-field  scanning;  optical 
communications;  optical  fiber;  optical  waveguide;  refracted  near-field  scan- 
ning;  refracted-ray   scanning;    resolution  limit. 


1.   INTRODUCTION 

Roughly  a   half-dozen  methods   have   been  proposed  and   implemented  for  determining   the 
refractive-index  profile  of  an  optical-fiber  waveguide.     At  least  half  of   these  will    no 
doubt  find  use   in  one  or  another  application. 

One  of  the  more  promising  techniques  for  a  simple  high-resolution  measurement  is  known 
as   refracted  near-field  scanning  or,   preferably,    refracted-ray   scanning.     First  proposed  and 
demonstrated  by   Stewart  [1],   the  method  was   further  developed  by  White   [2].      In  addition,    I 
have  analyzed  the  precision  of  the  method  [3,4]. 

The  purpose  of  this  Technical    Note  is  to  describe  refracted-ray  scanning  in  detail,   to 
analyze  it  as  a  measurement  system,   to  document  a  particular  system  for  implementing  it,   and 
to  show  some  exemplary  results.     First,   however,    let  us  examine   some  of   the  alternate 
approaches   that  are  either  in  use  or  have  been  proposed. 

Certain  methods   illuminate  the  fiber  in  a  direction   perpendicular  to  the  axis;   that  is, 
they  look  through  the  fiber  sides  rather  than  at  its  end  face.     These  transverse-illumina- 
tion methods  have  the  advantage  that  they  can  be  used  to  examine  a  fiber  in  real   time,   as  it 
is  being  produced.      Likewise,    these  methods  may  be  used  to  determine  the  fiber's  properties 
at  several    points  without  having  to  cut  the  fiber  to  prepare  an  end. 

One  of   the   simplest  techniques   is   that  developed  by   Liu,   who  used  an  oil-immersion 
method  similar  to  that  used  in  crystallography   [5].     By   placing  a  short  piece  of  fiber  in  a 
microscope  and   varying   the   index  of  a  matching  fluid,    Liu  was  able   to  measure  the   indices  of 
the  core  and  cladding  of  a  step  fiber. 

*Electromagnetic  Technology  Division,   National    Engineering  Laboratory 


Later  workers  illuminated  immersed  fibers  transversely  with  a  laser  and  placed  the 
apparatus  in  an  interference  microscope  or  interferometer  [6-11].  The  index  difference 
between  core  and  cladding  can  be  obtained  fairly  readily  and  without  detailed  calculation 
[6-8],  as  can  the  alpha  parameter  of  a  fiber  that  is  assumed  to  have  a  power-law  profile 
[9].  In  addition,  formulas  can  be  derived  for  calculating  the  full  index  profile  from 
measurements  of  the  fringe  intensity  [10,11].  The  overall  accuracy  of  these  measurements  is 
on  the  order  of  10  percent,  and  spatial  resolution  is  not  great  because  the  finite  width  of 
the  fringes  limits  resolution. 

Transverse-scattering  methods,  including  backscatter,  have  been  implemented  by  several 
groups  [12-14].  With  these  techniques,  a  fiber  (also  Immersed  in  oil)  is  illuminated  trans- 
versely by  a  laser,  and  the  far-field  diffraction  pattern  is  observed.  A  great  deal  of  data 
is  required,  so  that  an  automatic  data-acquisition  scheme  is  almost  mandatory;  the  index 
profile  is  calculated  by  computer.  The  accuracy  of  the  method  is  greater  than  of  transverse 
interferometry,  but  no  simple  or  direct  methods  for  determining  alpha  or  core-cladding  index 
difference  have  been  developed. 

More  recently,  Marcuse  and  Presby  have  developed  a  focusing  method  that  uses  transverse 
illumination  and  may  be  used  to  examine  either  fibers  or  preforms  [15-17].  This  technique 
scans  the  near  field  of  the  transversely  illuminated  fiber  and  calculates  the  index  profile 
from  data  acquired  by  a  computer.  Because  circular  symmetry  is  assumed  in  the  calculation, 
the  resulting  profile  is  necessarily  symmetric,  even  if  the  fiber  is  not. 

Other  methods  illuminate  the  fiber  along  the  axis;  either  a  thin,  specially  prepared 
sample  or  a  flat,  properly  prepared  end  is  required.  While  not  amenable  to  real-time  mea- 
surement as  a  fiber  is  drawn  nor  to  nondestructive  testing,  these  methods  provide  a  more 
direct  evaluation  of  the  index  profile,  and  certain  of  them  offer  resolution  that  is  com- 
parable with  conventional  optical  microscopy. 

Axial  or  slab  interference  microscopy  (as  opposed  to  transverse  interference  micro- 
scopy) requires  a  thin  slice  of  fiber  to  be  placed  in  an  interference  microscope  or  inter- 
ferometer [18,20].  The  interferogram  is  generally  interpreted  by  computer.  In  any  case, 
the  slab  is  treated  purely  as  a  phase  object  for  the  purpose  of  interpreting  the  inter- 
ference pattern;  because  a  graded-index  fiber  has  focusing  properties,  it  may  not  be  regard- 
ed as  a  phase  object  unless  the  sample  is  extremely  thin,  say  of  the  order  of  tens  of  micro- 
meters. Further,  it  is  difficult  to  determine  what  errors  arise  from  the  effects  of  polish- 
ing; these  may  include  changes  of  the  surface  composition  of  the  material  and  deviation  from 
flatness  as  a  result  of  the  polishing  rate's  varying  with  glass  composition.  As  with  trans- 
verse interferometry,  the  number  of  fringes  is  small,  so  the  spatial  resolution  is  apt  to  be 
relatively  low. 

Fresnel  reflection  from  the  cleaved  or  polished  end  of  a  fiber  has  been  used  to  measure 
the  index  profile  [21,23].  This  method  is  direct  and  precise  and  has  high  resolution.  It 
is  difficult  because  the  changes  in  reflectance  are  typically  less  than  1  percent;  highly 
stable  radiometric  techniques  must  be  employed.  This  method  has  recently  been  improved  with 
an  index-matching  technique  [22].  However,  the  reflectance  of  the  face  may  change  as  the 
result  of  contamination  or  exposure  to  the  atmosphere. 


other  methods  that  have  been   discussed  in  the  literature   depend  on  examination  of   the 
far-field  interference  pattern   [24-26],    on   special    properties   of   self-focusing  fibers   [27] 
and  on  the  spatial    coherence  of  the  light  emerging  from  an  optical    fiber  [28].      (See  also 
Ref.    29  for  a  more  detailed  review  of   index-profile  measurements.) 

2.      NEAR-FIELD   SCANNING 

Near-field  scanning  depends  on  the  existence  of  a  local    acceptance  angle   (or  numerical 
aperture)   at  any  point  on  the  entrance  face  of  the  fiber  [30].     That  is,    the  vertex  angle  of 
the  cone  that  can  enter  the  fiber  and  be   guided  depends  on   the   index  of   the  core  at  the 
point  of   illumination. 

The  fact  was   first  exploited  by   Sladen,   Payne,   and  Adams   [31,32].     This   group   illumi- 
nated the  end  of  a   short  fiber  with  a  uniform,   lambertian   source.     They   scanned  the  exit 
face  of  the  fiber  with  a  detector  and  thereby  generated  a  plot  of  the  output  power  as  a 
function  of  position  along  a   diameter  of   the  fiber. 

The  exact  proportionality  between  the  local    numerical   aperture  and  the  power  coupled 
into  the  fiber  holds  only   for  guided  rays.     Therefore,   Adams,   Payne,   and  Sladen  had  to 
derive  correction  factors  for  the  leaky-ray  contribution  to  the  power  that  propagates  into 
the  fiber  [33,34].     With  these  correction  factors,   they  are  able  to  calculate  the  index 
profile  to  a  high  degree  of  accuracy. 

Other  workers  devised  a  similar  scheme  for  measuring  the  index  profile  [35,36].     Rather 
than  illuminate  the  entire  entrance  face  of  the  fiber,   these  groups  focused  a  laser  beam  to 
a  point  and  observed  the  total   power  radiated  at  the  far  end  of  the  fiber  as  a  function  of 
the  position  of  the  point  on  the  entrance  face.     This  variation  of  near-field  scanning  also 
requires  correction  for  leaky  modes;    in  principle,    it  differs  only   slightly  from  the  ori- 
ginal   technique.* 

In  practice,  entrance-face  scanning  may  be  preferable  to  because  an  absolute  calibra- 
tion  is   possible.      In  addition,    imperfections  or  mode  coupling  along   the  fiber  will    have 
minimal    effect  on   the   relative  powers   transmitted  by   excitations  at  different  points  on   the 
entrance  face.     The  same  is  not  necessarily  true  of  the  exit-face  scan   (the  original   meth- 
od);  although   I  know  of   no  detailed  studies   to  this  effect,    it  seems  quite  possible  that 
such  factors  as  mode  coupling  will    distort  the   results   in   an  unpredictable  way,   particularly 
in  a  fiber  with  fine  structure  like  the  characteristic  index  dip  along  the  axis   (see  Section 
8,   below). 

There   is  another  difficulty  with  both  methods  of   near-field  scanning.     This   is   the 
problem  of  making  the  leaky-mode  corrections  with  confidence.     Microbending  or  slight  devi- 
ations of   the  fiber  from  circularity  may  drastically  alter  the  leaky-mode  propagation;    in- 
deed,  Petermann   has   suggested  that  the  leaky-mode  correction  may  be  unnecessary   for  fibers 


*For  simplicity  of   nomenclature  and  for  avoiding  confusion,    I   shall    refer  to  both   these 
methods  as  near-field  scanning.     When   it   is  necessary   to  distinguish  between   them,    I   shall 
use  the   terms  entrance-face   scanning  and  exit-face   scanning.      Likewise,    I  will    not  use  the 
term  refracted  near-field  scanning,    but  rather  will    refer  to  this  method  as  ref racted-ray 
scanning. 


with   nearly   square-law  profiles   [37].     Thus,   whereas  near-field  scanning   is  both   simple  and 
elegant,    it  may  be  that  its  major  use  will    be  to  give  excellent,  qualitative  information 
about  the   index  profile. 

3.      REFRACTED-RAY   SCANNING 

Stewart's  innovation  was  to  illuminate  a  fiber  with  a  focused  beam  whose  vertex  angle 
greatly  exceeds  the  acceptance  angle  of  the  fiber  and  to  observe  not  the  rays  that  are 
trapped  by  the  fiber,   but  the  rays  that  are  refracted  at  the  core-cladding  interface  [1]. 
With  this  method,   it  is  often  possible  to  eliminate  entirely  the  deleterious  effects  of  the 
leaky  modes  and  to  develop  a  system  that  generates  the  index  profile  directly,  without  the 
need  for  a  correction  factor. 

To  understand  the  principle  of  the  method,   it  is  helpful    to  examine  Fig.   1.     In  that 
figure,  a  fiber  whose  index  is  n  is  surrounded  by  a  liquid  whose  index  is  n^.     The  fiber  is 
illuminated  by  a  focused  cone  of  rays;   the  ray  that  is  drawn  in  the  figure  represents  the 
most-oblique  ray  in  the  cone,  a  marginal    ray  for  the  lens  that  focuses  the  beam.     Because 
the  angle  of  incidence  exceeds  the  acceptance  angle  of  the  fiber,   the  ray  is  refracted  both 
at  the  entrance  face  of  the  fiber  and  at  the  periphery. 
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Figure   1.     The  lower  the   index  of  the 
fiber,   the  greater  the  vertex  angle  of 
the  emergent  cone. 


Suppose  now  that  we  alter  the  index  of  the  fiber  in  some  way,   for  example  by  changing 
the  fiber  entirely.     The  ray  that  emerges  is  now  represented  by  the  dashed  line  in  Fig.   1; 
as  the  index  of  the  fiber  is  decreased,   the  angle  formed  by  the  marginal    ray's  intersection 
with  the  axis  increases  accordingly. 

Now  we   introduce  an  opaque,   circular  stop  directly  behind  the  fiber;   the  stop   is 
designed  to  intercept  all   but  the  outermost  rays  of  the  emergent  cone.     A  conical    shell    is 
transmitted  beyond  the  stop.      In  the  plane  of   the  stop,   the   inner  radius  of  the  shell    is 
always  equal    to   the   radius  of   the  stop.     The  outer  radius,   however,    varies  with  the   index  of 


the  fiber.  Therefore,  the  power  that  is  transmitted  around  the  stop  also  varies  with  index. 
The  heart  of  refracted-ray  scanning  lies  in  the  fact  that,  under  the  right  conditions,  the 
change  of  transmitted  power  is  precisely  proportional  to  the  change  of  index  of  the  fiber. 

Figure  1  used  a  uniform,  unclad  fiber  for  tutorial  purposes.  In  practice,  all  that  is 
required  is  a  fiber  whose  index  variation  is  a  function  of  radius  only.  The  vertex  of  the 
incident  cone  is  scanned  across  the  entrance  face  of  the  fiber;  the  power  transmitted  by  the 
stop  is  a  linear  function  of  the  index  of  the  fiber  at  the  point  of  illumination. 

4.   STEP  FIBER 

For  illustrative  purposes,  we  begin  by  examining  a  meridional  ray  in  a  step  fiber.  The 
extension  to  the  general  case  is  straightforward. 

Figure  2  traces  the  path  of  the  marginal  ray  that  is  focused  onto  the  face  of  the 
fiber.  The  fiber  is  immersed  in  a  fluid  whose  index  is  n,  .  The  index  of  the  fiber  core  is 
n;  that  of  the  cladding  is  n  . 


Figure  2.  The  path  of  a  meridional  ray 
through  a  step  fiber. 


The  ray  enters  the  system  at  angle  o  to  the  axis  of  the  fiber.   It  suffers  a  couple  of 
refractions  through  a  microscope  coverslip  (whose  thickness  is  not  shown)  before  reaching 
point  2.  Because  the  coverslip  makes  a  right  angle  to  the  axis  of  the  system,  we  may  easily 
show  that 


sin  e 


n.  sin  9„. 


;i) 


n|_  sin    e^  is   the   numerical    aperture   (NA)    of   the  cone   that  converges  onto  the  face  of   the 
fiber  and   is  equal    to  the  NA  of   the  microscope  objective  used  to  focus  the  light. 
If  we  apply  Snell's   law  at  points  2  and  3,   we   find  that 


and 


n.  sin  e^  =  n  sin  qL   >  (2) 


n  sin  e,  =  n  cos  e^  .  (3) 


If  we  square  and  add  these  equations  and  use  Eq.  (1),  we  find  that 


sin^  ^  "^  "c  si"^  ^3  ^  "^  •  ^^^ 


We  may  eliminate  the  term  in   9'   by  applying  Snell's  law  at  points  4  and  5: 


n     sin  e:.  =  n,    sin   e. 


and 

sin  9'   =  n,    cos   6^  .  (6) 

Finally,  we  square  and  add  Eqs.    (5)   and  (6)   and  combine  with  Eq.    (4)   to  learn  that 

n^  =  n^  +  sin^   9  -  sin^   9' .  (7) 

Equation   (7)   shows  how  the  vertex  angle   9'   of  the  emergent  conical    shell   of  rays  is 
related  to  the  core  index  n  of  the  fiber. 

5.     A  MORE-GENERAL  CASE 

Consider  the  general   case  of  a  skew  ray  incident  on  a  graded-index  fiber  with  angle  of 
incidence  9  and  arbitrary  azimuthal   angle  <i>;   this  is  shown  in  Fig.   3,  with  the  cover  slip 
omitted  for  clarity. 

Snell's  law  may  be  written  as 

n  sin  i  =  constant,  (8) 

where  i   is  the  angle  of  incidence  at  an  arbitrary  interface.     Equation   (8)   is  valid  whether 
the  index  change  is  sudden  or  gradual.     If  X  is  the  vacuum  wavelength  of  light,   then  the 
magnitude  k  of  the  wavevector  in  a  medium  is  Z-nn/X;   therefore,   Snell's  law  may  also  be 
written   in   the  equivalent  form, 

k   sin  i   =  k     =  constant.  (8a) 


Figure  3.  The  path  of  a  skew  ray  through  a  graded- index  fiber.  6  is  the  angle  of  inci- 
dence, and  (f),  which  lies  in  a  plane  perpendicular  to  the  axis,  is  the  azimuthal 
angle.  The  coverslip  at  the  front  of  the  cell  is  suppressed  for  clarity. 


That  is,  the  component  of  the  wave  vector  perpendicular  to  the  index  gradient  is  conserved 
as  the  ray  propagates.  White  has  used  this  principle  to  derive  a  general  form  of  Eq.  (7). 
It  is  applicable  to  any  fiber  whose  index  variations  are  purely  radial. 

Let  us  number  points  1  through  5  as  in  Fig.  2,  remembering,  however,  that  the  ray  path 
inside  the  fiber  is  in  general  curved  and  need  not  lie  in  a  plane.  Using  Eq.  (8a),  we  find 
that 


sin  9  =  n.  sin  G„  =  n(r)  sin  ei 


:9) 


where  n(r)  is  the  index  of  the  fiber  at  a  distance  r  from  the  axis  of  symmetry.  From  point 
2  until  point  4,  the  index  variations  become  purely  radial,  so  the  axial  component  of  the 
wave  vector  remains  constant;  that  is. 


n(r)  cos  e„'  =  n,  cos  6i- 

c  L  0 


(10) 


Finally,   at  point  5, 


n.    sin   Br  =  sin 


(11) 


Combining  Eqs.  (9)  through  (11)  as  before,  we  find  that 


?  7  7  7 

n  (r)  =  n.   +  sin  9  -  sin  e' 


(12) 


This  is  the  same  as  Eq.  (2)  of  Ref.  2.  When  n(r)  is  nearly  equal  to  n^  we  may  rewrite  Eq. 
(12)  in  the  form 

2n  An(r)  =  sin  6  -  sin  6'.  (12a) 

An(r)  =  n(r)  -  n,  «  n(r),  n ,;  it  is  not  the  delta  parameter  of  the  fiber.  We  discuss  the 
importance  of  Eqs.  (12)  and  (12a)  in  the  following  sections. 

6.   RADIOMETRIC  ANALYSIS 

A  lambertian  source  is  one  whose  radiance  is  independent  of  angle.  Many  thermal 
sources  approximate  lambertian  sources;  lasers  and  many  light-emitting  diodes  do  not.  We 
begin  this  section  by  assuming  that  the  source  is  lambertian  and  later  generalize  to  a  non- 
lambertian  case. 

Into  a  cone  whose  vertex  angle  is  e' ,  a  lambertian  source  emits  total  power  proportion- 
al to  sin^  6'  [38].  In  the  ref racted-ray  technique,  we  insert  an  opaque  stop  behind  the 
fiber  to  block  some  of  the  rays  refracted  by  the  fiber,  as  in  Fig.  1.  If  the  source  is  lam- 
bertian, then  the  power  transmitted  around  the  stop  may  be  expressed  as 

P(6')  =  A  (sin^  e'  -  sin^  6  ),  ,  (13) 

where  6^  is  the  angle  subtended  by  the  stop  and  A  is  a  constant  of  proportionality. 

If  we  let  n(r)  =  n^_  (no  fiber  in  the  system),  then  we  find  from  Eq.  (12)  that  9'  =  6. 
If  we  call  the  power  transmitted  around  the  stop  Pq,  then  we  find  that 

A  =  P^/(sin^  9  -  sin^  e^).  (14) 

Finally,  if  we  combine  Eqs.  (12a),  (13)  and  (14),  we  find  that 

P( 9' )  -  P 
2  n^  An(r)  =  - — p    °  (sin^  9  -  sin^  9^).  (15) 

0 


All  the  terms  in  this  equation  are  constants  except  P(9')  and  An(r).  Thus,  when  An  is 
small,  the  power  that  propagates  beyond  the  stop  is  in  principle  precisely  proportional  to 
the  index  profile  of  the  fiber.   (The  derivation  has  been  made  for  guided  rays  only;  see 
Ref.  2  and  below  for  a  discussion  of  leaky  rays.) 

Equation  (15)  could  be  used  for  calibration  of  the  system  if  the  two  angles,  9  and  9^, 
could  be  measured  with  sufficient  precision.  In  particular,  it  is  difficult  to  measure  9^, 
because  this  is  the  angle  subtended  by  the  stop  at  a  point  inside  the  fiber-- that  is,  at  the 
virtual  image  of  the  point  source  as  seen  from  the  rear  of  the  cell.   In  part  for  this  rea- 
son, it  may  be  preferable  to  devise  a  direct  calibration  method  (see  Section  11,  below). 


Many  sources  are  not  lambertian.  For  example,  a  laser  beam  that  is  expanded  and  passed 
through  a  small -diameter  lens  more  nearly  approximates  a  uniform  point  source  than  a  lamber- 
tian source.  In  contrast,  edge-emitting  diodes  and  semiconductor  lasers  emit  their  radia- 
tion more  strongly  in  the  forward  direction  than  does  a  lambertian  source.  Such  sources  can 
be  approximated  at  least  roughly  by  writing  their  radiance  L(9)  in  the  form, 

L(e)  =  Lq  cos*"  0  ,  (16) 

where  we  assume  circular  symmetry.  A  lambertian  source  is  described  by  m  =  0  and  a  uniform 
point  by  m  =  -1.  Other  sources,  such  as  edge-emitting  diodes,  may  be  approximated  by  using 
Eq.  (16)  with  m  an  integer  between,  say,  5  or  10,  depending  on  the  specific  source. 

As  noted  above,  the  power  radiated  by  a  lambert  source  into  a  cone  is  proportional  to 
sin^e;  see,  for  example,  Eq.  (3.17)  of  Ref.  35.  If  we  follow  the  derivation  leading  to  that 
equation,  we  find  that  the  fraction  of  the  total  power  radiated  into  a  cone  is 

1  -  cos"^^^  8,  (17) 

which  reduces  to  sin^  9  when  m  =  0.  (This  relationship  and  all  subsequent  ones  are  also 
valid  when  m  =  -1.)  Equation  (13)  becomes,  in  the  general  case, 

P(e')  =  A  (cos"^^^  0^  -  cos'"^^  6')  ,  (18) 

and  Eq.  (14)  becomes 

A  =  Pq/(cos'^"'^  e^  -  cos"'''^  6)  .  (19) 

If  Eq.  (18)  is  to  describe  the  fiber  index  profile,  it  must  be  linear  in  sin^  9', 
because  Eq.  (12a)  shows  that  An(r)  is  linear  with  sin^  6'.  I  have  not  found  it  fruitful  to 
try  to  relate  Eq.  (18)  directly  with  sin^  8;  however,  it  is  instructive  to  plot 
(1  -  cos"^^  8)  [Eq.  (17)]  as  a  function  of  sin^  8.  This  is  shown  in  Fig.  4  for  several 
values  of  m. 

The  most-important  cases  are  the  lambertian  source  and  the  uniform  point  source,  for 
which  m  is  equal  to  0  and  -1.  The  former  case  is  obviously  linear,  and  the  latter  is  nearly 
so  for  numerical  apertures  up  to  0.6  or  more.  Fortunately,  what  is  important  is  that  the 
curves  be  linear  only  over  the  small  range  of  angles  defined  by  the  stop  diameter  and  the 
numerical  aperture  of  the  focusing  lens. 

As  m  increases,  the  curves  become  less  linear.  Also,  they  deviate  substantially  from 
the  line  m  =  0.  Therefore,  unless  m  is  precisely  0,  a  calibration  based  on  Eq.  (15)  will  be 
inval  id. 
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Figure  4.     A  function   useful    for 
describing  the  power  transmitted  by  a 
non-1 ambertian   source. 


7.      LEAKY-RAY   ANALYSIS 

Suppose  that  a  circular  fiber  is  illuminated  off  axis  by  a  focused  cone  of  rays.     Geo- 
metric optics  predicts   that  certain   rays  will    be  trapped  by   the  fiber,   whereas  certain 
others  will   be  refracted  out  of  the  fiber;  we  refer  to  these  cases  as  guided  and  refracted 
rays   (or  bound  and  radiation  modes). 

Snyder  and  his  colleagues  have   shown   that  certain   rays  are  not  truly  guided  but  never- 
theless suffer  high  loss  and  therefore  are  greatly  attenuated  after  transmission  along  a 
relatively   short  length  of  fiber;   they  called  these  leaky   rays   [39-41].     The  presence  of 
leaky   rays   invalidates  a   guided-ray  analysis   such  as  that  of   Sections  4,   5  and  6. 

Adams,   Payne  and  Sladen  have  derived  correction  factors  to  be  used  with  the  near-field 
scanning  technique   [33,34].      However,    as   noted  above,   there   is   some  doubt  as   to  the  utility 
of  these  correction  factors   in   some  cases.     The   ref racted-ray  technique  in  many  cases  avoids 
the  necessity   for  a  correction  factor  by   adjusting  the  angular  subtense  of  the   stop   so  that 
the  leaky-ray  power  transmitted  by   the  fiber  is   intercepted  by   the  stop;   only   truly 
refracted  rays  are  passed  by  the  stop.     When   the  leaky-ray  contribution   is  excluded,   the 
guided-ray  analysis  becomes   valid   [33]. 

If  a  fiber  is  excited  at  a  point  at  a  distance  r  from  the  center,   leaky  rays  will   be 
excited  only   for  angles  of   incidence   9  given   by   [2,41] 


(n^(r) 


n^^)   /   (1 


2  2  2 

p     cos     (J))    >  sin 


>   n^ir)   -  n^2, 


(20) 
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<|)  is  the  azimuthal  angle  shown  in  Fig.  3,  a  is  the  core  radius  and  p  is  (r/a);  all  other 
terms  have  been  defined  before. 

Following  White,  we  exclude  leaky  rays  by  restricting  angles  of  incidence  to  values 
larger  than  the  largest  value  of  the  left  side  of  Eq.  (20).  That  is,  we  make  the  opaque 
stop  large  enough  that  rays  pass  the  stop  only  if  their  angle  of  incidence  is  larger  than 
that  given  by 

sin^  9  >  (n^(r)  -  n^2)/(l  -  p^).  (21) 

The  right  side  of  this  equation  defines  the  largest  angle  of  incidence  at  which  leaky  rays 
will  be  excited  at  a  given  radius. 

If  we  apply  Eq.  (11)  to  Eq.  (21),  we  find  that  leaky  rays  will  be  excluded  provided 
that  6  exceeds  the  value  of  e'  given  by 

sin^  6^  >  sin^  6'  =  n^^  -  n^(r)  +  [(n2(r)  -  n^)/[l  -   p^)]  .         (22) 

White  has  applied  Eq.    (22)   to  a  fiber  with  a  power-law  profile  and  quoted  his   results 
[2].     A  power-law  profile  may  be  described  by  the  equations 

n^(r)   =  n^^   (1   -  2p"a)    ,  (23a) 

and 

n^^  -  n^(a)   =  n^^   (1   -  2a)    ,  (23b) 

where 

n^  =   n(0)   and  A  =    (n^^  -  n^^)/2nj-.  (23c) 

If  we  use  these  equations  in  Eq.  (22),  we  find  after  a  few  lines  of  algebra  that 

sin^  e^  >  (n^^  -  n^^)  +  2nQ2  p^A  (1  -  p°')/(l  -  p^),  (24) 

which   is   the   same  as  White's  Eq.    (10)    [2]. 

To  eliminate  the  leaky   rays  entirely,   we  must  examine  Eq.    (24)   when  the  p  term  on   the 
right  side   is  largest.     Physically,    it  is  clear  that  this   term  must  be  a  maximum  when   p  =   1, 
because  leaky   rays  are  excited   in  greater  numbers  as  we  move  away  from  the  axis   toward  the 
core-cladding   interface.     We  may  also   show  this  analytically   by   letting  R  =  p"^  and  differ- 
entiating F(R) ,   where 


i 


F(R)  =  R(l  -  R^)/(l  -  R)  (25) 

and  x  =  a/2.  The  derivative  is  equal  to 
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dF/dR  =  [1   -   (l+x)    R^  +  X  r'^"^^]!!   -  R)'^  (26) 

which   is   indeterminate  when  R  =   1.     We  may  evaluate   the  limit  by   substituting  R  =   1   -   e  and 
letting  e  ^  0.      We  find  in  this  way  that  the  derivative   is   indeed  0  when  R  =  1,   and  that 
F(R)    is  a  maximum  at  the  core-cladding  boundary. 

Returning  now  to  Eq.    (24),   we  evaluate  F(R)    in   the   limit  R  ^  1   by  applying  I'HOpital's 
rule.     The  calculation  is  straightforward  and  leads  to  the  conclusion  that  F(R)   approaches 
a/2  as  R  approaches   1. 

Thus,   the  angular  subtense  of  the  stop  must  exceed 

sin^   6     >    (n^^   -n^)+n^aA,  (27) 

which  is  equivalent  to  Eq.  (13)  of  Ref.  2. 

If  the  parameters  in  Eq.  (27)  are  reasonably  well  known,  then  it  is  an  easy  matter  to 
calculate  the  minimum  stop  diameter.  In  the  event  that  all  the  parameters  are  not  known, 
then  it  is  possible  to  measure  the  first  term  on  the  right  side  of  Eq.  (27).  Let  I  be  the 
axial  acceptance  angle  of  the  fiber.  Then,  according  to  Ref.  2  or  Ref.  38,  Section  2.8, 

sin^  I  =  n„^  -  x\}   .  (28) 

0     c 

We  now  return  to  Eq.  (11)  and  replace  G  and  e'  with  I  and  I',  where  I'  is  the  angle  at  which 
the  guided  ray  emerges  from  the  cell.  Combining  Eqs.  (11)  and  (27)  in  this  way  shows  that 

sin^  r  =  n,  ^  -  n  ^  .  (29) 

L     C 

2     2  7 

Thus,  the  term  (n,   -  n  )  in  Eq.  (27)  is  equal  to  sin*^  I',  where  I'  is  the  inner  vertex 

angle  of  the  hollow  cone  of  light  that  is  refracted  by  the  fiber.  An  alternate  way  of 

writing  Eq.  (27)  is  therefore  to  replace  (n[_^  -  n^^)  with  sin  I'. 

Finally,  ^.n^^A   is  equal  to  sin^  I,  where  sin  I  is  sometimes  called  the  theoretical 

numerical  aperture  of  the  fiber.  Therefore,  Eq.  (27)  may  be  rewritten 

sin^  e^  >  sin^  I'  +  (a/2)  sin^  I.  (30) 

The  stop  must  subtend  an  angle  larger  than  the  value  of  e'  given  by  Eq.  (27)  or  (30). 
By  way  of  example,  let  us  choose  a  fiber  for  which  n  =  1.45  (approximately  the  value  for 
vitreous  silica),  sin  I  =  0.25,  and  a  =  2.3.  Assume  further  that  n[_  is  one-half  percent 
larger  than  n^,  (it  is  preferable  to  choose  n|_  >  n^  to  avoid  exciting  cladding  modes).  I 
have  chosen  these  values  deliberately  to  overestimate  sin  9. 

Using  Eqs.  (29)  and  (30),  we  find  that  sin^  I'  =  0.021  and  therefore  that  sin  6^  >  0.3. 
A  microscope  objective  with  a  numerical  aperture  in  the  neighborhood  of  0.5  will  be  requir- 
ed; fortunately  these  are  quite  common. 
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A  profile  with  an  index  step  may  be  thought  of  as  a  power-law  profile  for  which  a 

approaches  infinity.  According  to  Eq.  (30),  there  is  no  value  of  6  for  which  leaky  rays 
can  be  completely  excluded  if  a  is  unbounded. 

When  r  =  0,  e'  is  a  minimum  (we  exclude  from  this  discussion  the  possibility  of  an 

index  dip  at  the  center).  Therefore,  to  allow  observation  of  the  complete  index  profile, 

the  angular  subtense  of  the  opaque  stop  must  not  exceed  the  value  of  e'  given  from  Eq.  (11) 
by 

7  7  7  7  7 

sin  6  <  sin  0'  =  sin  e  +  n,   -  n  ,  (31) 

s  L     0 

which   is  just  Eq.    (11)   with   r  =  0.     We  may  write   n     in  terms  of  A  by  manipulating  Eq.    (23c), 
which  yields 

n/  -  n/  .  Np2,  (32) 

where  Np  =  sin  I  =  2n^A.  The  second  equality  holds  because  a  is  very  small.  Using  Eq. 
(32)  in  Eq.  (31),  we  find  that 

sin^  e  <  sin^  6  +  n,  ^  -  n  ^  -  N^^  ,  (33) 

s  L     C     F 

which  is  White's  Eq.  (15)  [2]. 

Equation  (31)  or  (33)  delimits  the  angular  Subtense  of  the  largest  stop  that  may  be 
used  with  a  given  fiber.  Equation  (24)  delimits  the  greatest  angle  of  refraction  at  which 
leaky  modes  will  be  found.  This  angle  increases  with  radius.  Following  White,  let  us 
assume  that  the  subtense  of  the  stop  is  precisely  that  given  by  Eq.  (33);  in  fact,  it  must 
be  slightly  less  than  this  value.  If  we  scan  a  fiber,  leaky  modes  will  not  appear  beyond 
the  stop  until  the  radius  is  sufficiently  large  that  the  angle  e'  defined  by  Eq.  (24) 
exceeds  the  subtense  of  the  stop  as  given  in  Eq.  (33).  To  find  an  implicit  expression  for 
this  radius,  we  set  the  right  sides  of  these  two  equations  equal  and  find  that 

N^/Np^  ^  (1  _  p(a+2)/(^  _  p2j  (34) 

where   p  =   (r/a),   and  N  =  sin   9  is   the  numerical    aperture  of  the   incident  cone. 

We  calculate  the   limit  of  Eq.    (34)   as   p  approaches   1   and  solve  for  a.      Because  p  =   1  at 
the  core-cladding   interface,   this  procedure  yields   the  largest  value  of  a  for  which   leaky 
modes  will    not  be  a   factor.     The  calculation   is  carried  out  using   I'HOpital's   rule;   the 
result  is 

a  =   2[(N^/Np2)    -   1]    ,  (35) 
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which  is  Eq.  (18)  of  Ref.  2.  For  a  fiber  for  which  sin  I  is  0.25  illuminated  by  a  micro- 
scope objective  whose  NA  is  0.5,  a  =  6;  a  increases  to  10.5  if  sin  I  decreases  to  0.2. 
Thus,  leaky  rays  can  be  handled  for  nearly  all  practical  fibers  except  possibly  step  fibers. 

For  step  fibers,  we  examine  Eq.  (34)  in  the  limit  as  a  ^  «>.  Because  p  is  at  most  1, 
the  right  side  simply  approaches  1/(1  -  p^).  Therefore, 

r/a  <  (1  -  Np^/N^)  (36) 

for  leaky  rays  to  be  excluded  by  the  stop.     This  relation  is  precise  only  when  e.  has  been 
optimized  for  each  fiber.     For  the  examples  used  in  connection  with  Eq.    (35),  p  =  0.87  and 
0.92,   in  that  order.     Beyond  these  radii,   an  indeterminate  fraction  of  the  leaky-ray  power 
is  blocked  by  the  stop;   therefore,   it  is  not  possible  to  calculate  a  correction  factor. 
Thus,   the  refracted-ray  method  is  not  useful    for  scanning  step  fibers  beyond  about  nine 
tenths  of  their  radius. 

8.      RESOLUTION   LIMIT,    EDGE   RESPONSE  AND  DEPTH  OF  FOCUS 

White  has  argued  that  the  spatial    resolution  limit  of  the  system  is  limited  by  the 
opaque  central    stop  [2].     The  stop  lies  roughly  in  the  entrance  pupil   of  the  condensing 
lens,  and  it  is  natural    to  project  the  stop  into  the  exit  pupil   of  the  microscope  objective, 
where  it  obscures  the  appropriate  fraction  of  the  exit-pupil    diameter.     Because  the  micro- 
scope is  diffraction  limited,   the  resolution  limit  may  be  expected  to  be  that  of  a  diffrac- 
tion-limited annular  aperture. 

This  argument  would  be  precisely  correct  if  the  condensing  lens  were  itself  a  diffrac- 
tion-limited imaging  lens;   in  that  case,   the  resolution  limit  in  the  image  plane  of  the 
condensing  lens  would  be  determined  by  the  annular  stop,  which  could  be  projected  to  any 
convenient  point  of  the  optical    system. 

Born  and  Wolf  have  calculated  the  Fraunhofer-dif fraction  pattern  of  a  circular  aper- 
ture, a  fraction  e  of  whose  diameter  is  obscured  by  an  opaque  circular  stop  [42];   it  is 

2 

Ke)    =    [2J^(B)/B  -    2e^J^(ee)/(e3)]^/(l-e^)      ,  (37) 

where  J^  is  the  Bessel   function  of  order  1,    B  =   i-n/x)   D  sin  e  is  the  normalized  dimension  in 
the  image  plane,  D  is  the  exit-pupil    diameter,   e  is  the  angle  by  which  the  observation  is 
off  axis,  and  the  intensity  at  the  center  of  the  pattern  is  1.     The  radius  of  the  Airy  disk 
(first  zero  of  the  Bessel   function)   of  a  clear  aperture  occurs  where  B  =  1.22Tr. 

Equation   (37)   describes  the  impulse  response  or  point-spread  function  of  a  lens  with  an 
annular  aperture.     When  we  scan  a  fiber,   we  translate  the  point-spread  function  at  right 
angles  to  the  index  variations  of  interest.     Therefore,  we  require  knowledge  of  the  edge 
response  of  the  lens,   rather  than  the  impulse  response  [43].     This  is  the  convolution  of  the 
impulse  response  with  a  unit  step  in  the  direction  normal    to  the  direction  of  translation; 
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Location  of  edge 


Figure  5.     Geometry  used  to  calculate  the  edge  response  of  a  diffraction-limited  lens  with  a 
central    stop. 


It  may  be  calculated  by  carrying  out  the  volume  integral    indicated  in  Fig.   5,   in  which  the 
sombrero-like  function  represents  Eq.    (37).     The  integration. 
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has  been  carried  out  numerically  for  values  of  z  between  0  and  0.9  [44].         is  the  angle 
shown   in  Fig.    5.     Although   the  edge   response   is   defined  for  an  opaque  edge,    it  should 
approximately  describe  the  response  of  an  optical    system  to  the  sharp  index  step  of  interest 
to  us. 

We  may  define  the  edge-response  width  in  a  manner  analogous  to  the  definition  of  rise 
time  or  transition  duration  in  electrical   engineering:     The  edge-response  width  is  the  dis- 
tance between  the  two  points  where  the  intensity  of  the  image  is  10  percent  and  90  percent 
of   the  maximum.     This  parameter  gives  a  measure  of   the   resolution   limit  of   the   system  as  a 
function  of  the  normalized  stop  diameter  e. 

The  solid  curve  of  Fig.    6  shows  the  edge-response  width  as  a  function  of  relative   stop 
diameter  e.     The  ordinate  is  normalized  to  the  Airy-disk  radius,   so  the  graph  has  universal 
applicability.     When  e  =  0,   the  edge-response  width  is  very  close  to  the  Rayleigh  limit; 
this   is   in   good  agreement  with  measurements   of   resolution  of  bar  targets.      (Incidentally, 
the  curve  lies  very  close  to  the  function  l/(l-e),   which   is  what   I   guessed  by   replacing  D 
with  D(1-g)   in  the  expression  for  the  Airy-disk  radius.) 
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Figure  6.     Resolution  limit  as  a  function  of  central    stop  diamter  for  the  case  of  a 
diffraction-limited  lens.     Solid  line  is  the  calculated  edge  response: 
experimental   points  are  measured  index-step  widths. 
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Figure  7.     Measured  edge  response  of  refracted-ray  scan  for  various  values  of  e. 
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To  measure  the  edge-response  width  of  the  ref racted-ray  scanning  system  (to  be  describ- 
ed shortly),  I  cleaved  some  home-made  quartz  fibers,  immersed  them  in  oil,  and  scanned 
across  the  edge  along  a  diameter  that  was  roughly  at  right  angles  to  the  direction  along 
which  the  fiber  had  been  cleaved.  Figure  7  shows  a  typical  set  of  such  scans  for  several 
values  of  e.  This  parameter  varies  by  about  10  percent  depending  upon  whether  the  light  is 
focused  into  the  fiber  or  the  oil;  the  numbers  in  Fig.  7  are  the  averages  of  these  two 
values. 

Data  derived  from  the  scans  of  Fig.  7  and  other  scans  are  shown  as  discrete  points  in 
Fig.  6.  Whereas  there  is  slight  effect  from  increasing  e,  the  data  do  not  agree  well  with 
the  theory  and  show  that  resolution  is  diffraction  limited  as  long  as  e  is  less  than  0.8. 
(These  data  at  least  roughly  agree  with  White's;  however,  White's  approximate  calculation 
does  not  agree  with  the  curve  of  Fig.  3,  which  was  generated  without  recourse  to  a  gaussian 
approximation  to  the  sombrero  function,  2Jj^(x)/x.) 

The  condensing  lens  is  merely  a  collector  of  light  and  does  not  project  a  high-quality 
image.  Indeed,  it  may  be  thought  of  as  a  part  of  the  detector,  for  the  system  would  be  un- 
changed if  the  condensing  lens  and  detector  were  replaced  by  a  single,  large  detector. 
Thus,  we  may  regard  the  refracted-ray  apparatus  as  a  scanning  microscope  that  uses  an  annu- 
lar detector. 

Kermisch  has  pointed  out  that  a  scanning  microscope  is  precisely  equivalent  to  a  con- 
ventional imaging  system,  provided  that  the  detector  be  replaced  with  an  incoherent  source, 
and  the  source  with  a  point  detector  [45].  Therefore  our  system  is  equivalent  to  a  micro- 
scope with  a  central  stop  in  the  condenser  optics;  its  resolution  limit  (and  other  charac- 
teristics) may  be  analyzed  by  performing  calculations  on  this  equivalent  system.  Considera- 
tions based  on  the  Abbe  theory  show  that  the  resolution  limit  ought  to  be  roughly  constant 
as  the  size  of  the  central  stop  is  increased.  (An  extreme  case  is  dark-field  microscopy; 
similar  considerations  show  that  resolution  is  comparable  to  that  of  conventional,  bright- 
field  microscopy  and  can  even  be  slightly  better  [46].)  At  any  rate,  the  resolution  limit 
of  conventional  microscopy  changes  by  no  more  than  30  percent  or  so  as  varying  conditions  of 
illumination  result  in  different  degrees  of  coherence  [42]. 

Change  of  the  stop's  angular  subtense  is  equivalent  to  a  similar  change  of  illumination 
in  conventional  microscopy.  Therefore,  we  expect  the  resolution  limit  of  our  system  to  vary 
only  slightly  with  changing  stop  position,  and  this  indeed  seems  to  be  the  case.  I  cannot 
at  present  account  for  the  sudden  loss  of  resolution  at  large  values  of  e.* 

Depth  of  focus  is  related  to  resolution  limit  and  is  one  factor  that  may  determine  the 
precision  with  which  the  end  of  the  fiber  must  be  perpendicular  to  the  axis  of  the  system. 
For  a  diffraction-limited  optical  system  such  as  a  microscope  objective,  depth  of  focus  may 
be  written  [47] 

DF  =  nX/2(NA)^.  (39) 

*  Footnote  added   in  proof:     W.   J.    Stewart  has   shown   that  the  resolution  limit  is  actually  a 

minimum  for  certain  values  of  e;    that   is,   resolution   limit  increases   as  c   is   decreased  below 

the  values  examined  in  this  paper.     Using  these   results,   Reid  and  Stewart  have   succeeded   in 

developing  a  system  with  0.35  \im  resolution  limit  [48,49]. 
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If  NA  =  0.5  and  A  -   633  nm,  the  depth  of  focus  DF  is  about  2  ym.  If  we  focus  on  the  center 
of  a  lOO-pm  diameter  fiber,  then  the  edges  will  be  in  sharp  focus  only  if  the  face  of  the 
fiber  deviates  from  perpendicularity  by  less  than  (2/50)  rad,  or  about  2°. 

9.   CONDENSING  LENS 


With  no  fiber  in  the  cell,  the  hollow  cone  transmitted  through  the  cell  will  have  a 
numerical  aperture  NA  equal  to  that  of  the  microscope  objective.  If  the  fiber  index  is 
slightly  less  than  that  of  the  matching  fluid,  then  the  cone  angle  is  increased  slightly  but  HI 
is  still  of  the  order  of  the  incident  cone  angle.  The  hollow  cone  must  be  focused  by  a  con- 
densing lens  onto  a  detector. 

Lenses  (other  than  microscope  objectives)  are  more  commonly  described  in  terms  of  their 
F  number  <))  than  numerical  aperture.  Figure  8  shows  the  geometry  used  for  calculating  the 
relationship  between  NA  and  F  number.  In  air,  the  NA  is  given  by 


NA  =  sin  U 


D/2(f'^  +  D^/4) 


V2 


(40) 


from  which  we  may  readily  calculate  that 

4.  =  f'/D  =  (1/2)(NA"^  -  D" 


h 


[41) 


Thus,  an  NA  of  0.5  is  equivalent  to  an  F  number  of  about  0.9.  (Table  1  compares  NA  and  F 
number.)  We  must  design  a  condensing  lens  with  an  effective  F  number  <^{1   +  magnification) 
less  than  1  [50]. 


Figure  8.  Vertex  angle  of  a  lens. 


Table  1 


F  Number,  Numerical  Aperture,  and  Vertex  Angle 


m_ 

0.15 
0.20 
0.25 
0.40 
0.50 
0.65 


3.3 

2.4 

1.9 

1.15 

0.87 

0.58 


~~ir 

8.6' 
12 
15 

23.5 
30 
41 
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Single-element  condensing  lenses  may  be  purchased  readily  with  F  numbers  of  2  or  some- 
what less.  To  achieve  an  effective  F  number  less  than  1,  we  must  design  a  compound  lens, 
the  first  element  of  which  is  operated  with  the  object  distance  about  one  half  the  focal 
length. 

To  minimize  the  aberrations  of  the  lens,  we  attempt  to  make  the  deviation  of  the  ray  by 
each  lens  a  constant;  likewise,  we  choose  the  elements  in  such  a  way  that  the  deviation  at 
each  surface  is  about  one  half  the  total  deviation  of  the  lens. 

Figure  9  shows  a  thin  lens  whose  focal  length  is  f  projecting  an  image  with  conjugates 
I   and  Ji' .     The  angles  u  and  u'  are  defined  as  shown,  as  is  the  deviation  6.  The  theorem 
concerning  the  external  angle  of  a  triangle  shows  that 


^^^^^^^ 

h 

1 
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1 
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V 

Figure  9.  Deviation  and  other  para- 
meters of  the  condensing  lens. 


6  =  u  +  u' 


(42) 


This  result  is  true  for  thick  lenses  as  well  as  thin.  In  paraxial  approximation. 


u  =  h/£  ,  (43a) 


and 


u'  =  h/£'  .   (43b) 


Using  Eq.  (42)  and  the  lens  equation,  we  find  immediately  that 

6  =  h/f.   (44) 

Now  consider  a  three-element  condensing  lens  (Fig.  10).  For  the  moment,  consider  the 
lenses  to  be  thin  and  in  contact;  all  have  focal  length  f.   If  we  demand  a  magnification 
m  =  2,  then  the  image  distance  l\   is  equal  to  7.%^^   where  i^   ^^  ^^^   object  distance.  If  the 
deviation  by  each  lens  is  6,  then  the  total  deviation  is  36;  therefore,  looking  at  the  lens 
as  a  whole,  we  find  that 
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36  -  (h/z^)   +  (h/2Jii) 


(45) 


Because  6  =  h/f  the  object  distance  £i  is  f /2, 
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Figure  10.  Conjugates  of  the  three 
element  condensing  lens. 


A  simple,  paraxial  ray  trace  using  the  lens  equation  shows  that  the  first  lens  projects 
a  virtual  object  to  the  primary  focal  point  of  the  second  lens;  the  second  lens  projects  the 
image  to  =»;  and  the  third  lens  projects  the  image  to  its  secondary  focal  point.  Thus,  to 
minimize  spherical  aberration,  the  lenses  should  be  planoconvex,  each  oriented  with  its 
curved  side  toward  its  own  long  conjugate.  (The  first  element  should  in  principle  be  menis- 
cus, but  such  a  lens  may  not  be  available.)  If  the  lenses  have  been  molded  and  are  to  be 
placed  into  a  barrel,  they  must  be  centered  and  have  their  edges  ground  by  an  optics  shop  to 
keep  other  aberrations  to  a  minimum.  Image  quality  is  still  poor,  as  with  all  condensing 
lenses,  so  a  uniform,  large-area  detector  is  required. 

In  reality,  the  elements  are  comparatively  thick  with  respect  to  their  focal  lengths. 
Therefore,  it  may  be  necessary  to  follow  the  path  of  a  marginal  ray  through  the  lens  to  be 
certain  that  one  of  the  central  elements  does  not  act  as  the  aperture  stop  for  the  system. 
Apart  from  this  caveat,  the  design  based  on  thin-lens  elements  is  probably  adequate.  The 
effective  focal  length  of  the  real  lens  will  exceed  the  calculated  value  somewhat. 

For  most  cases,  a  three-element  lens  working  at  the  magnification  of  2  will  probably  be 
adequate;  if  not,  a  four-element  lens  may  be  necessary.  It  is  easy  to  generalize  Eq.  (44) 
to  the  case  of  an  N-element  lens  with  magnification  m.  The  result  of  an  analysis  exactly 
similar  to  the  preceding  is 


(1  +  1/m)  f'/N  , 


(46) 


£|;j  =  (1  +  m)  f'/N  , 


(47) 


where  a'    is  the  image  distance. 

Finally,  the  design  of  the  condenser  is  complicated  and  could  lead  either  to  a  very 
large  lens  or  to  a  thick  lens  with  barely  sufficient  clearance  between  the  lens  and  the  end 
of  the  capillary.  W.  C.  Meixner  of  Valtec  suggested  the  use  of  plastic  Fresnel  lenses 
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instead.  Aspheric  Fresnel  lenses  are  available  from  more  than  one  manufacturer  with  F 
numbers  as  low  as  0.7.  These  lenses  are  corrected  for  spherical  aberration,  generally  pro- 
vided that  one  conjugate  is  at  infinity.  Therefore,  two  lenses  with  equal  diameter  (but  not 
necessarily  equal  focal  length)  will  be  required:  The  first  roughly  collimates  the  beam, 
and  the  second  focuses  onto  the  detector.  The  F  number  of  0.7  implies  that  the  numerical 
aperture  of  the  microscope  objective  may  exceed  0.5  if  necessary.  Much  of  the  experimental 
work,  including  the  final  calibration  run  and  fiber  scans,  reported  in  the  remainder  of  this 
paper  was  carried  out  after  a  Fresnel-lens  condenser  had  been  installed. 


10.  APPARATUS 

The  experimental  apparatus  is  sketched  in  Fig.  11.  The  HeNe-laser  beam  is  expanded  by 
the  first  microscope  objective  MO  and  focused  by  the  second  40X  microscope  objective  through 
a  microscope  cover  slip  onto  the  fiber.  The  fiber  is  held  in  the  moveable  cell,  whose  cross 
section  is  shown  in  Fig.  12  and  which  is  discussed  below.  The  fiber  passes  through  a  hole 
in  the  opaque  stop  OS.  The  refracted  rays  escape  the  fiber  and  are  focused  by  the  condens- 
ing lens  CL  onto  a  uniform,  1-cm-diameter  silicon  detector.  The  lens  is  designed  as 
described  in  the  previous  section;  the  elements  have  127-mm  focal  length  and  9-cm  diameter 
and  are  fixed  in  an  aluminum  housing  as  shown  in  Fig.  13.  Because  the  laser  is  polarized, 
the  x/4  plate  is  used  to  reduce  the  dependence  of  various  reflectances  on  angle. 

Figure  14  shows  two  rear  views  of  the  apparatus,  with  the  lens  and  detector  removed. 
The  cell  is  constructed  of  aluminum.  The  glass  window  is  fastened  with  a  good  grade  of 
epoxy.  The  cover  slip  is  fastened  to  a  mounting  plate  with  a  fast-setting  epoxy  that  may  be 
removed  with  a  razor  blade  if  the  cover  slip  has  to  be  replaced;  the  plate  is  screwed  to  the 
cell  housing  in  case  the  cover  slip  must  be  removed  and  cleaned.  The  fiber  is  introduced 
through  a  0.25-mm  capillary  tube  that  is  fastened  to  the  window  with  the  same  rapidly  set- 
ting epoxy,  so  that  the  capillary  tube  may  be  removed  easily  with  a  heat  gun.  The  outer  end 
of  the  tube  has  been  formed  into  a  funnel  to  guide  the  fiber.  The  stop  is  made  of  aluminum 
by  turning  it  on  a  lathe  and  is  fixed  in  place  with  piano  wire,  as  shown  in  Fig.  15.  Stray 
light  that  passes  through  the  hole  in  the  stop  is  intercepted  by  a  piece  of  black  tape  fast- 
ened to  the  center  of  the  lens. 

In  addition,  I  illuminate  the  fiber  with  a  white  light  placed  directly  behind  the  cell 
(or  sometimes  at  the  far  end)  and  inspect  the  entrance  face  with  the  beam  splitter  and  eye- 
piece assembly  shown  in  Fig.  11.  The  green  filter  is  used  to  eliminate  stray  laser  light. 
Once  the  beam  splitter  and  eyepiece  are  fixed  in  place,  they  may  be  used  to  locate  the  fiber 
and  focus  on  it  by  manipulating  the  cell  position.  This  procedure  is  accurate  to  better 
than  25  m;    I  adjust  the  focus  further  by  scanning  across  the  index  dip  or  other  sharp  fea- 
ture and  examining  the  profile  for  sharpness.  The  white-light  source  at  the  far  end  of  the 
fiber  may  be  replaced  with  a  detector  for  transmitted  near-field  scanning. 

The  prototype  system  used  an  ordinary  three-axis  translator  to  hold  the  cell.   I  scan- 
ned horizontally  by  hand  and  attached  the  micrometer  to  a  potentiometer  by  means  of  an  0 
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Figure   11.     Experimental    setup.     Q--quarter-wave  plate. 


--microscope  objective.      BS--beam 


splitter.     A--aperture  stop.     OS--opaque  stop.     F--fiber.     WL--white-light 
source.     CL--condensing  lens.     D--detector.     EP--eyepiece.     GF--green  filter. 
MY--experimenter. 


microscope 
objective 


window 


fiber 


Figure  12.     Detail   of  cell   assembly  and  focusing  optics.     Outside  the  capillary,   the  fiber 
is  confined  to  the   plane  perpendicular  to  the  direction  of  scan. 
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Figure  13.  Photograph  of  condensing  lens. 


i 
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Figure  14a. 


Figure  14b. 


Figure  14.  Photographs  of  the  cell  and  translation  stage,  condensing  lens  removed. 
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Figure  15.     Photographs  of  the  opaque  stop. 


ring  used  as  a   drive  belt;    the  potentiometer   is  used   in  a  voltage-divider  circuit  with   two 
batteries  or  a  plus-and-minus  15-V  power  supply  and  drives  the  horizontal   axis  of  the 
recorder.     This  system  is  easy  to  set  up  and  may  be  adequate  for  many  purposes.     Because  the 
resolution  of   the  micrometers   is  25   pm   (0.001   in),   the  fiber  may  be  positioned  to  a  preci- 
sion of  about  6   pm;   reproducibility   is   somewhat  poorer  than  this   figure. 

I   improved  the  prototype  system  by  employing  a  1-rpm  ac  motor  to  drive  the  horizontal 
scan  of  the  three-axis  manipulator.     With  a  toothed  wheel    and  matching  drive  belts,   this 
system  was  sufficiently  precise  to  measure  the  edge-response  curves  discussed  above.     For 
this  purpose,   I   replaced  the  wire-wound  10-turn  potentiometer  with  a  continuous  1-turn 
potentiometer. 

The   system   I  currently  employ  uses  a  three-axis  manipulator  whose  precision   is  2.5   um 
(0.0001   in)   or  better.     I  adjust  the  vertical   and  focusing  axes  by  hand  and  drive  the  hori- 
zontal   scan  with  the  clock  motor  and  toothed  drive  belts.     In  addition,   I  have  made  provi- 
sion to  scan  with  a  translation  stage  driven  by  a  digital    stepping  motor  with  0.2  pm  steps. 
Such  a  digital    drive  may  be  useful    for  computer  analysis,   but  the  resolution  limit  of  0.2   m 
is  insufficiently  fine  for  much  of  the  work  reported  here. 

With  one  exception,   the  optical   alignment  of  the  apparatus  is  not  critical.     I  intro- 
duce the  laser  beam  by  means  of  two  mirrors   (not  shown).     After  rough  alignment,   I  center 
the  opaque  stop  by  eye;   the  aperture  stop  creates  pronounced  diffraction  rings  that  make  the 
visual   alignment  fairly  precise.     After  introducing  the  condensing  lens  and  positioning  the 
detector  by  maximizing  its  output,    I  use  the  mirrors   to  maximize  the  power  through  the  sys- 
tem;  this  ensures  that  the  laser  beam  is  centered  about  the  aperture  stop  that  precedes  the 
microscope  objective.     The  opaque  stop  may  be  adjusted  further  at  this  time,   but  such  ad- 
justment does  not  seem  necessary. 

The  location  of  the  aperture  stop,   however,    is  critical.     The  aperture   stop  determines 
the  axis   of  the  cone  that  focuses  onto  the  face  of  the  fiber.      If   that  axis   is   not  very 
nearly  parallel    to  the  axis  of  the  fiber,   then  the  scans  will    be   skewed  slightly.      In  a 
graded-index  fiber,   the  cladding  index  shows  a   slight  upward  or  downward  slope  that  is  con- 
sistently  in  the  same  direction  on  both  sides  of   the  core.     To  correct  for  this  error,    I 
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found  that  the  aperture  stop  must  be  positioned  horizontally  to  an  accuracy  of  25  m   or 
so.  With  a  40X  microscope  objective  (4-nm  focal  length),  this  corresponds  to  a  pointing 
accuracy  of  about  5  mrad  or  one  third  of  a  degree.  The  adjustment  is  best  made  by  scanning 
a  fiber  repeatedly  until  the  cladding  or  the  Oil  levels  on  either  side  of  the  scan  are 
equalized.   In  addition,  to  avoid  creating  an  angular  error  in  the  horizontal  plane,  I 
introduce  the  fiber  by  holding  it  in  place  directly  above  the  capillary,  rather  than  to  one 
side. 

11.  CALIBRATION 

Data  are  taken  on  an  xy  recorder.     The  output  of  the  photodetector  passes  through  a 
trans-impedance  amplifier   (current-to-voltage  converter)   whose  output  enters   the  vertical 
amplifier  of   the   recorder.      I   checked  the  overall    linearity  of   the  electronic   system  by 
first  calibrating  a  set  of   neutral -density   filters  at  633  nm  on  a  spectrophotometer. 
Introducing   the  filters   into   the  beam  ahead  of  the  first  MO  yielded  the  curve   shown   in 
Fig.    16.     Electrical    linearity   is   therefore   good  over  a  factor  of  at  least  5. 
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Figure  16.  Electrical  linearity  of  the  system,  determined  by  placing  calibrated  filters 
into  the  laser  beam. 


White  [2]  has  calibrated  his  system  by  translating  the  stop  to  different  positions 
along  the  axis.  I  sought  a  more-direct  calibration  scheme,  which  might  be  more  amenable  to 
analysis  [3].  Ideally,  one  or  more  fibers  with  known  core  and  cladding  indices  would  make 
calibration  an  easy  matter.  Unfortunately,  such  a  set  of  fibers  is  not  available;  neither 


25 


is  there  a  guarantee  that  the  index  of  a  fiber  is  precisely   that  of   the  preform  from  which 
the  fiber  was  drawn. 

Malitson  has  shown  that  vitreous  silica  is  manufactured  with  sufficient  purity  that  the 
index  variation  from  sample  to  sample  or  manufacturer  to  manufacturer  is,   for  our  present 
purposes,   negligible  [51].     He  has  also  derived  Sellmeier  coefficients  that  allow 
calculation  of   the   index  at  any  wavelength  up  to  3.71   pm.     However,   Fleming  has  criticized 
the  use  of  Malitson's  numbers  for  calculating  the  index  of  a  fiber;   fibers  are  not  annealed, 
but  rather  are  chilled  or  quenched  and  may  not  have  the  same  index  as  the  annealed,   bulk 
material    [52,53].     Hence,   Fleming  has  determined  the  Sellmeier  coefficients  relevant  to 
quenched  silica. 

At  the  wavelength  632.8  nm,   the  Malitson  formula  gives  the  result  that  the  index  of  the 
silica  at  20°C   is   1.457018;   the  Fleming  formula  gives   1.457334  at  23.5°C   (see  Table  2). 
Malitson  quotes  a  temperature  coefficient  of  10.0  x  10"°  K"-*-;   thus  at  23.5°C,  the  index  of 
vitreous  silica  is  1.457053.     At  this  temperature,   the  Malitson  and  Fleming  values  differ  by 
0.00028,  an  amount  that  is  significant  when  compared  with  the  accuracy  to  which  the  indices 
of  the  immersion  fluids  are  known.     I  chose  the  average  value  of  1.45719   (not  1.45726  as 
incorrectly  reported  in  Refs.   3  and  4)   and  assigned  an  error  of  ±0.00014. 


Table  2 


Sellmeier  coefficients  for  vitreous  silica  (0.21-3.71  urn  and  20°C)  [49]  and  for  quenched 
silica  (0.44-1.53  m   and  23.5°C)  [50].  All  wavelengths  are  expressed  in  micrometers. 


fused  silica 
quenched  silica 


0.6961663 
0.696750 


0.4079426 
0.408218 


0.8974794 
0.890815 


fused  silica 
quenched  silica 


0.0684043 
0.069066 


0.1162414 
0.115662 


9.896161 
9.900559 


A^X 


2   2 

X  -A, 


2        2 
2  "^  72 


x^-x. 


X  -X, 


For  calibration,   I   used  several   fibers  drawn  by  hand  from  vitreous-silica  rods  as  well 
as  the  naked  core  of  a  plastic-clad  silica  fiber.      (Although  the  latter  was  more  convenient, 
I  wanted  to  verify  that  it  was  truly  silica.)     I  compared  these  fibers  with  each  of  five 
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Figure  17.  Calibration  scans  with  unclad  PCS  fiber.  Structure  near  right  edge  is  a  result 
of  the  fracture  and  does  not  represent  true  index  data. 


oils  sequentially.  The  indices  of  the  oils  were  supplied  by  the  manufacturer,  who  claimed 
an  accuracy  of  ±0.0005  [54]. 

Figure  17  shows  a  typical  calibration  run  with  the  silica  core.  When  the  laser  is 
focused  into  the  oil  rather  than  the  fiber,  the  rays  should  emerge  parallel  to  their  origi- 
nal direction;  we  might  expect  the  power  transmitted  to  be  independent  of  the  index  of  the 
fluid.  The  slight  variation  in  Fig.  17  may  be  accounted  for  by  the  change  in  reflectance  at 
various  oil -glass  interfaces.  Because  of  this  variation,  it  is  necessary  to  normalize  the 
results  to  the  power  transmitted  through  the  oil  alone. 
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Figure  18.  Normalized  voltage  as  a  function  of  index  difference  between  silica  and  oil,  for 
a  particular  stop  position.  Estimates  of  random  and  systematic  error  ^.r^ 
indicated. 
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Figure  18  shows  the  results  of  some  calibration  runs  similar  to  those  reported  in  Ref. 
3.  The  vertical  axis  is  (V^  -  Vq)/Vq,  where  V  is  the  voltage  read  by  the  xy  recorder  and  s 
and  0  stand  for  silica  and  oil.  The  horizontal  axis  is  n  -  n  ,  where  n  is  index  of 
refraction. 

The  greatest  barrier  to  a  linear  calibration  may  be  ensuring  the  absence  of  vignetting 
when  the  highest-index  oil  is  in  the  cell  (and  the  emergent  cone  was  largest).  It  is 
important  to  verify  the  linearity  of  the  system  over  a  large  range  of  indices  in  case,  for 
example,  some  of  the  fibers  to  be  tested  have  cladding  indices  less  than  that  of  vitreous 
silica.  Because  of  the  possibility  of  vignetting,  or  other  effect,  an  extrapolation  may  be 
dangerous  if  linearity  over  a  sufficient  range  has  not  been  demonstrated. 

A  calculation  along  the  lines  described  by  Nlatrella  showed  linearity  at  the  99  percent 
confidence  level  [55].  However,  I  did  not  calculate  a  conventional  line  of  best  fit  because 
the  points  have  uncertainty  in  the  horizontal  direction  as  well  as  the  vertical.  Rather,  I 
used  a  method  similar  to  another  suggested  by  Natrella.  The  center  of  mass  of  all  the 
points  as  well  as  the  means  of  each  of  the  four  sets  are  calculated.  Three  line  segments 
join  adjacent  means,  and  their  slopes  are  calculated.  The  best-fit  line  is  taken  to  be  the 
one  that  passes  through  the  center  of  mass  and  whose  slope  is  equal  to  the  average  of  the 
slopes  of  the  three  line  segments.  This  is  the  line  shown  in  Fig.  18. 

The  vertical  error  bar  in  Fig.  18  is  an  estimate  of  the  instrumental  limit  of  error  (or 
resolution)  of  the  xy  recorder;  the  actual  vertical  errors  are  also  of  this  order.  The 
horizontal  error  is  the  sum  of  the  index  uncertainties  just  mentioned,  as  well  as  an  addi- 
tional ±  0.0001  owing  to  the  effect  on  the  oils  of  the  uncertainty  of  the  oil  temperature 
measurement.  Because  the  large  horizontal  error  bar  is  not  reflected  in  large  vertical 
scatter,  I  conclude  that  most  if  not  all  of  the  horizontal  error  (except  the  temperature 
component)  is  a  systematic  error  that  is  common  to  each  oil  and  has  approximately  constant 
magnitude  and  sign;  the  two  horizontal  error  bars  show  the  random  and  systematic  errors 
separately.  (See  Table  3  for  a  compilation.) 

When  running  an  unknown  sample,  I  usually  choose  the  oil  whose  index  at  25°C  is  about 
1.464;  generally  this  value  is  slightly  higher  than  the  cladding  index,  and  the  values  of 
AV/V  for  the  scan  remain  within  the  range  defined  by  Fig.  18.  In  the  unusual  event  that 
some  of  the  values  fell  outside  that  range,  the  oil  would  have  to  be  changed  to  avoid  having 
to  extrapolate  the  line  in  Fig.  18. 

After  running  the  scan,  I  remove  the  fiber  to  determine  the  voltage  corresponding  to 
the  index  of  the  oil  without  the  complicating  effect  of  the  shadow  of  the  fiber.  (It  is 
equivalent  to  continue  to  scan  for  several  fiber  radii  beyond  the  cladding.)  I  measure  the 
voltage  difference  corresponding  to  the  core-cladding  index  difference  and  compute  aV/V  for 
the  fiber.  The  value  can  be  converted  directly  to  index  difference  by  using  the  slope  of 
the  calibration  curve  of  Fig.  18.  For  relative  measurements,  it  is  not  necessary  to  know 
the  precise  index  of  the  oil;  therefore.  An  should  be  accurate  to  +0.0006,  as  detailed  in 
Table  3.  Because  it  eliminates  the  index  of  the  oil  from  the  measurement,  this  procedure  is 
more  accurate  than  calibrating  each  measurement  with  a  vitreous  silica  fiber  and  monitoring 
the  temperature  of  the  oil. 
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The  approximate  index  of  the  fiber  is  needed  to  calculate  the  numerical  aperture, 
2 
NA  -   /(2n  An).  However,  because  n  is  large,  an  error  in  n  is  not  nearly  so  severe  as  an 

error  in  An,  and  the  NA  may  be  calculated  precisely  with  only  approximate  knowledge  of  the 

index  of  the  fluid. 


Table  3 


Error         Budget 


Vertical 


Instrumental 
limit  of  error 


±0.01 


random 


Temperature  coefficient 

of  refractive  index,  +0.25  K 


Horizontal 
+0.0001 


random 


Index  of  silica 


Index  of  oils  at  25°C 


±0.00014 
±0.0005 


systematic 

assumed 
systematic 


SUBTOTALS 


±0.0001 
±0.00064 


random 
systematic 


Projection  of 
vertical   error 


±0.0005 


random 


TOTALS 


±0.0006 
±0.0006 


random 
systematic 


The   index  of  any  of  the  oils   is   significantly  temperature  dependent;    for  the  D  line, 
the  coefficient  is  of   the  order  of  4  x  10"^  K".     Thus,   with   slight  increase   in  complexity, 
we  could  calibrate   the   system  using  only  a  single  oil    and  varying   its  temperature  by  a  few 
tens  of  kelvins  while  measuring  to  a  precision  of  a  few  tenths  of  a  kelvin. 

With  this  scheme,  we  could  get  many  more  points  along  the  calibration  curve,  but  would 
have  to  measure  the  index  of   the  oil    as  a  function  of   temperature  at  precisely   the  wave- 
length of  interest.     Further  thoughts  on  index  measurements  of  liquids  are  detailed  in 
Appendices  A  and  B. 
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12.      MEASUREMENTS  ON  ACTUAL  FIBERS 

In  this  section,  I  report  the  results  of  refracted- ray  scans  of  a  number  of  actual 
communication  fibers,  on  a  comparison  with  conventional  near-field  scans,  and  on  a  com- 
parison with  another  laboratory. 

I  obtained  a  step  fiber,  designated  DF-1 ,   that  had  been  subjected  to  an  exit-face  scan 
and  a  far-field  scan  by  other  workers  in  this  laboratory  [56].     I  performed  additionally  a 
refracted-ray  scan  and  an  entrance-face  scan  on  this  fiber. 

Figure  19a  shows  the  refracted-ray  scan.     The  scan  reveals  a  central    structure  that  is 
actually  an  index  dip  surrounded  by  a  nearly  circular  index  dip   (as  may  be  verified  with  the 
eyepiece).     In  addition,   the  scan  shows  that  the  fiber  has  a  low-index  cladding  layer 
surrounded  by  what  is  evidently  silica. 

Incidentally,   if  we  knew  that  the  outermost  layer  of  a  fiber  were  in  reality  silica, 
then  we  could  use  the  index  of  that  layer  as  a  calibration  level.     Most  of  the  fibers  I  have 
tested  had  cladding  indices  that  were   (within  experimental   accuracy)   equal   to  that  of  sili- 
ca;  some,   however,   had  substantially  lower  index.     It  is  therefore  unwise  to  presume  any- 
thing about  the  fiber  in  the  absence  of  a  priori   knowledge  that  it  is  clad  with  pure  silica. 

The  apparently  high-index  region  near  the  edge  of  the  core  of  the  fiber  DF-1  is  an 
artifact  that  results  from  the  presence  of  leaky  rays.     That  the  structure  is  not  real   may 
be  verified  by  translating  the  opaque  stop  either  in  or  out  and  noticing  that  the  magnitude 
of  the  peaks  varies  with  stop  position   (not  shown  in  a  figure).     These  artifacts  may  be 
eliminated  by  translating  the  stop  axially  to  the  proper  position,   but  at  the  expense  of 
having  to  recalibrate  the  system  and  probably  losing  resolution. 

We  may  also  use  the  data  to  calculate  the  core-cladding  index  difference  An  and  employ 
Eq.    (36)   to  determine  the  portion  of  the  core  that  is  free  of  leaky  rays;   this  region  in 
indicated  in  Fig.   19a  and  corresponds  well   with  the  appearance  of  the  artifacts. 

Figure  19b  shows  exit-face  scans  of  both  a  short  length  and  a  kilometer  length  of  this 
fiber,   both  performed  at  857   nm.     The  upper  scan  nearly   resolves  the  index  dip  but  suffers 
from  the  leaky-ray  problem.     The  lower  scan,   on  the  other  hand,   is  free  of  leaky  rays  and 
suggests  that  the  fiber  core  does  indeed  have  a  nearly  constant  index  of  refraction.     How- 
ever,  the  index  dip  has  been  greatly  exaggerated  (compared  with  Fig.   19a),   probably  because 
of  mode  coupling  by  scattering  or  diffraction  out  of  the  core  center  as  a  result  of  the  long 
propagation  distance. 

Figure  19c   shows  entrance-face  scans  of  short  and  long  pieces  of  this  fiber.     To  make 
these  scans,   I   simply  removed  the  condensing  lens  from  the  system  and  placed  the  output  end 
of  the  fiber  into  an  oil-filled  capillary  tube  and  brought  it  nearly  into  contact  with  the 
detector.     As  before,   the  scan  of  the   short  length  shows   the  index  dip  with  what  is  evi- 
dently its  proper  magnitude  but  suffers  from  the  leaky-ray  problem  near  the  core  edge.     The 
scan  of  the  longer  piece,   however,   distorts  the  index  dip  substantially,   not  only  by  exag- 
gerating its  magnitude,   but  also  by   injecting  a  spike   into  the  center  of  the  pattern.      (This 
spike  is  easily  explained  by  noticing  that,   because  of  the  oscillatory  nature  of  the  index 
dip,   the  center  of  the  fiber  contains  a  tiny  annular  waveguide.     Apparently  this  waveguide 
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Figure  19a.     Refracted- ray  scan  of  step  fiber.     Arrowheads  show  region  calculated  to  be  free 
of  leaky  rays. 


Figure   19b.      Exit-face   scans  of   step  fiber  at  857   nm.      Upper--short  length;   lower--kilometer 
length. 
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Figure  19c.     Entrance-face  scan  of  step 
fiber  at  633  nm.     Upper--short  length; 
lower- -kilometer  length. 


has  lower  loss  than  the  fiber  as  a  whole.     The  scan  does  not  resolve  the  waveguide,  which 
appears  as  a  sharp  spike.) 

The   refracted-ray   scan   shows  details  of  the   step  fiber  that  are   not  shown  or  are   dis- 
torted by  conventional    near-field  scanning.      However,    the   refracted-ray   scan   suffers  from 
the  presence  of  leaky  rays  near  the  periphery  of  the  step-fiber  core.     It  is  therefore 
prudent,  when  examining  step  fibers,   to  supplement  the  scan  with,   say,   entrance  face  scans 
taken  with  different  numerical   apertures. 

To  assess  the  leaky-ray  problem  in  graded-index  fibers,    I  performed  a  refracted-ray 
scan  and  an  entrance-face   scan  of  a  56-cm  length  of   fiber  DF-D;    these  are  shown   in  Fig. 
20.      (Because  the   refracted-ray   scan  has  to  be   inverted  to  make  the  positive  n  axis  verti- 
cal,   the  upper  scan   is   in  the  opposite  direction  from  the  lower  two.)     The  horizontal    scales 
are  the  same  for  the  two   scans;   the  vertical    scales   have  been  adjusted  to  give  nearly  the 
same  amplitude.     The  cladding  of  fiber  DF-D   is  not  uniform,   as   I   have  determined  by  scanning 
several    freshly  cleaved  ends;   the  structure   in  the  cladding  of  Fig.    20   (upper)    is   real    and 
not  an  artifact  of  the  apparatus.     Figure  20   (lower)    is  an  exit-face  scan   of  a   long  piece  of 
the  same  fiber;   the  horizontal    scale  is  not  the  same  as   those  of  Figs.    20   (upper)   and 
(center). 
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Figure  20.     Scans  of  a  graded-index  fiber.     Llpper--refracted-ray  scan;  center--entrance-face 
scan,   short  length;   lower--entrance-face  scan,  kilometer  length. 

Even  a  casual    inspection  shows  the  two  near-field  scans  to  be  substantially  more 
rounded  than   the   refracted-ray   scan. 

Among  the  other  fibers   I   scanned  were  graded-index  fibers  designated  DF-A,   B,   D  and 
E.     These  fibers  had  been  used  in  a  round  robin  to  intercompare  radiation-angle  or  numeri- 
cal-aperture measurements  made  in  each  of  several    laboratories  [57].     Figure  21    shows  a  scan 
of  one  of  these  fibers,   DF-B.     This  curve  is  chosen  deliberately  to  show  the  glitch  to  the 
left  of  the  center,  which  is  probably  the  result  of  a  speck  of  dirt  or  a  shard  of  glass  too 
small    to  be  seen  with  the  eyepiece. 


Figure  21.     Refracted-ray   scan  of  fiber  DF-E.     Glitch  to  left  of  center  is  likely   due  to 
contamination  by  dirt  or  shard  of  glass. 
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Figure  22.  Comparison  of  numerical - 
aperture  measurements  by  round  robin 
(RR)  with  calculation  based  on  index- 
profile  measurements  (NBS). 
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Figure  23.  Core-cladding  index  differ- 
ence measured  by  NBS  and  by  other  labo- 
ratory (OL)  that  also  uses  refracted- 
ray  scanning.  Average  difference 
between  measurements  is  3  percent. 
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The  maximum  acceptance  angle  U  for  meridional   rays  is  given  implicitly  by  the  formula 

(48) 


2  2  2 

sin^  U  =  n„^   -  n/ 
0  c 


sin  U  is  sometimes  called  the  numerical  aperture  of  the  fiber.  Measured  n^  and  n^  and 
calculated  sin  U  from  Eq.  (48).  Table  4  compares  the  calculated  values  with  the  means  of 
the  values  obtained  for  each  fiber  by  the  laboratories  in  the  round  robin.  (The  same  data 
are  shown  graphically  in  Fig.  22.) 

Table  4.  Comparison  with  a  Round  Robin 


Numerical  Aperture 


Fiber 

DF-A 

DF-B 

DF-D 

DF-E 


Round-robin  average  [56] 
0.235  ±  2.9%* 
0.163  +  2.4% 
0.202  ±  1.6% 
0.195  ±  0.9% 


Calculated,   this  report 
0.264 
0.174 
0.208 
0.202 


^One  standard  deviation,  expressed  in  percent  of  mean. 


In  all   four  cases,   the  values  calculated  from  index  data  are  slightly  larger  than  the 
sample-average  numerical   aperture.     The  participants  in  the  round  robin  were  instructed  to 
measure  the  angle  between  the  points  at  which  the  intensity  was  5  percent  of  the  peak, 
whereas  the  calculated  numerical   aperture  will    include  all   the  light  transmitted  by  the 
fiber.     Therefore,  we  expect  the  calculated  values  to  exceed  slightly  the  measured  values, 
and  this  is  precisely  the  case. 

Finally,   I  was  fortunate  to  be  able  to  make  an  informal   comparison  with  M.  J.   Saunders 
of  Bell    Laboratories,  who  also  uses  refracted-ray  scanning  [58].     Owing  to  the  round  robin, 
which  was  going  on  at  the  same  time,  we  were  easily  able  to  exchange  a  total   of  six 
fibers.     We  measured  the  index  differences  shown  in  Table  5  and  Fig.   23.     The  other  labora- 
tory's value  is  generally  slightly  higher  than  that  of  NBS,   but  the  arithmetic-average 
difference  between  our  measurements  is  3  percent  even  though  we  used  different  calibration 
techniques. 
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Table  5.  Interlaboratory  comparison  of  index  difference  An. 


Fiber 

DF-A 
DF-B 
DF-D 
DF-E 
MJ-A 
MJ-B 


This  report 
0.0237 
0.0104 
0.0148 
0.0139 
0.0150 
0.0156 


Other  laboratory 
0.235 

o.ono 

0.0156 
0.0140 
0.0154 
0.0163 


Average 


Difference 
-0.8X 
+5.8 

5.4 

0.7 

2.7 

4.5 
+3% 


13.   ADDITIONAL  REMARKS 

One  of  the  major  practical  problems  with  refracted-ray  scanning  is  the  accumulation  of 
dirt  particles  or  shards  of  glass  on  the  inside  of  the  microscope  cover  slip.  Probably 
these  particles  are  deposited  by  the  fiber  onto  the  inside  of  the  capillary;  they  are  later 
picked  up  by  other  fibers  and  deposited  in  turn  onto  the  cover  slip.  In  any  case,  such  par- 
ticles accumulate  on  the  cover  slip  and  eventually  interfere  with  the  measurement. 

The  problem  may  be  less  severe  in  a  laboratory  environment,  where  relatively  few  fibers 
are  tested  than,  say,  in  a  production  environment.  I  simply  clean  the  fibers  in  acetone 
after  cleaving  them;  when  the  cover  slip  becomes  contaminated  with  particles,  I  wash  or 
replace  it. 

K.  I.  White  of  the  British  Post  Office  informs  me  that  they  insert  their  fibers  by 
using  a  tube  within  a  tube,  rather  than  a  single  capillary.  They  push  the  fiber  through  a 
fine  tube,  cleave  it,  and  place  the  tube  and  the  fiber  in  an  ultrasonic  cleaner  (which  con- 
tains the  same  oil  as  the  cell).  They  then  withdraw  the  fiber  into  the  tube  and  insert  the 
tube  into  a  larger  tube  that  is  fixed  to  the  cell.  They  feel  that  much  of  the  contamination 
problem  is  eliminated  in  this  way. 

White  further  informed  me  that  their  cell  is  cylindrical  and  mounted  in  a  V  block  for 
easy  removal  and  accurate  replacement.  Such  a  design  will  greatly  simplify  the  removal  and 
cleaning  of  the  cover  slip. 

Another  problem  is  that  the  fiber  drifts  slightly  within  the  cell.  The  rate  of  drift 
varies,  but  is  often  as  large  as  a  micrometer  or  so  per  minute.  With  multimode  fibers,  this 
is  not  important  unless  we  are  interested  in  examining  the  index  dip  or  other  fine  features; 
however,  a  micrometer  of  drift  may  cause  significant  error  in  a  scan  of  a  single-mode  fiber. 

The  drift  can  be  arrested  by  clamping  the  fiber  just  outside  the  cone  of  refracted 
rays;  I  simply  use  an  alligator  clip,  suitably  padded  with  a  soft  rubber  gasket  (not  shown 
in  the  figures).  Because  the  total  motion  of  the  capillary  is  only  100  pm  or  so,  the 
alligator  clip  need  not  be  fixed  to  the  translation  stage.  With  the  fiber  thus  clamped, 
drift  is  very  slight,  whereas  without  clamping,  the  fiber  sometimes  drifts  unacceptably 
before  a  full  scan  can  be  completed. 
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APPENDIX  A.     MEASUREMENT  OF   INDEX  OF  REFRACTION 

It  may  be  desirable  to  measure  the  index  of  the  matching  fluids  directly,  for  example 
if  the  calibration  is  to  be  carried  out  by  varying  the  temperature  of  a  single  fluid.     If  a 
refractometer  is  not  available,   the  measurement  may  be  carried  out  using  one  of  several 
methods  that  employ  a  prism  spectrometer.     The  theory  and  practice  of  these  methods  are 
described  in  detail   by  Longhurst  [59].     In  this  appendix,   I  will   discuss  several   factors 
that  limit  the  precision  of  liquid-prism  and  total -internal -reflection  methods  of  measuring 
the  index  of  a  liquid.     I   include  this  material    here  because  I  have  been  unable  to  find 
similar  calculations  elsewhere. 

There  is  no  need  to  repeat  Longhurst' s  description  of  the  minimum-deviation  method. 
The  index  is  determined  from  the  formula, 

-  sin   (g  +  9)/2  ,^^. 

"  sin  a/ 2         '  ^'^^' 

where  a   is  the  prism  angle  and  e  is  the  angle  of  minimum  deviation.  At  minimum  devia- 
tion, a  +   6  =  21.,  where  l^   is  the  angle  of  incidence  at  the  first  prism  face,  and  a  ray 
travels  through  the  prism  perpendicular  to  the  bisector  of  angle  a.  To  analyze  the  errors 
inherent  in  this  method,  we  first  differentiate  with  respect  to  both  a  and  e.  The  results 
are 

An  =  (n/2)[cot  a/2  -  cot  (a+6/2]Aa  (A2) 

a 

and 

An^  =  (n/2)[cot  (a+6)/2]Ae,  (A3) 

where  An  is  the  error  that  results  in  n  if  a  is  known  to  an  accuracy  of  Aa  and  where  An.  is 
defined  similarly.  Aa  and  ao  are  instrumental  errors  and  may  be  taken  as  one  half  the 
finest  divisions  on  the  divided  circles  of  the  spectrometer  (see  also  Ref.  60). 

For  mathematical  convenience,  we  may  treat  these  errors  as  correlated  and  add  their 
magnitudes  to  estimate  the  total  instrumental  error.  Paying  careful  attention  to  sign,  we 
find  that  Eqs.  (A2)  and  (A3)  are  positive  as  written;  therefore,  the  sum  An^  of  the 

instrumental  errors  An  and  An.  is 

a       0 

An^  =  (n/2)  Aa  cot  (a/2),  (A4) 

where  we  assume  that  AS  and  Aa  are  equal.  (Near  minimum  deviation,  9  is  a  slowly  varying 
function  of  angle  of  incidence.  The  orientation  of  the  prism  need  be  accurate  no  better 
than  several  minutes  to  measure  e  precisely.) 
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For  a  l-minute  spectrometer,  Aa  =  0.5'  or  about  0.15  mrad;  if  a  =  60°  and  n  =  1.5,  An, 
is  consequently  about  0.0002.  This  is  about  the  precision  required  for  index-profile 
measurements. 

Taken  by  itself,  Eq.  (A4)  suggests  that  a  should  be  made  as  large  as  possible,  provided 
that  it  is  not  so  large  that  light  will  not  pass  through  the  prism.  Figure  Al,  however, 
shows  that  as  a  increases,  the  deviation  becomes  a  stronger  function  of  angle  of 
incidence.  As  a  result,  measurement  of  the  minimum  deviation  e  becomes  more  dependent  on 
the  precise  orientation  of  the  prism.  Inasmuch  as  increasing  a  from  60°  to  80°  can  be 
expected  to  reduce  experimental  error  by  only  about  30  percent,  little  or  nothing  will  be 
gained  by  making  a  much  larger  than  the  common  value  of  60°. 
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Figure  Al.  Deviation  as  a  function 
90     of  angle  of  incidence  for  various 
prism  angles. 


The  prism  must  also  be  levelled;  that  is,  the  edge  of  the  prism  must  be  made  parallel 
to  the  axis  of  rotation  of  the  prism  table.  For  the  degree  of  accuracy  required  here, 
levelling  error  is  potentially  significant,  as  may  be  shown  by  referring  to  Fig.  A2.   In 
that  figure,  the  yz  plane  is  the  plane  that  contains  the  face  of  the  prism.  The  incident 
ray  OB  strikes  the  prism  with  angle  of  incidence  i  =  AOB.  However,  because  OB  is  inclined 
to  the  xy  plane  with  angle  e  =  COB,  the  spectrometer  will  measure  angle  i'  =  AOC.  We  wish 
to  calculate  the  difference  between  the  measured  value  i'  and  the  true  value  i. 

We  construct  a  vector  OB  whose  magnitude  is  1  and  apply  the  law  of  cosines  to  the 
triangle  OAB: 


I 


2  2    ? 

X  +  1  -  2x  cos  1  =  x  +  y 


(A5) 
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prism  face 


Figure  A2.  Geometry  for  calculation  of  error  resulting  from  improper  leveling  of  prism. 

X,  y,  and  z  are  the  line  segments  indicated  on  the  figure  and  y^  +  z^  is  the  length  of  AB. 
We  express  x  and  y  in  terms  of  the  angles  e  and  i',  and  then  demand  that  e  be  sufficiently 
small  that  cos  e  =  1  and  sin  e  =  e.  Eq.  (A5)  therefore  becomes 


2  2      2 

cos  i'  +1-2  cos  i'  cos  i  =  sin  i'  +  e 


(A6) 


When  e  =  0,  Eq.  (A6)  may  be  written  in  the  form, 


2  2 

cos  i  +  1  -  2  cos  i  cos  i  =  sin  i  . 


:a7] 


We  may  combine  Eqs.    (A6)   and   (A7)   by   subtracting  and  writing  the   result  in  the  form. 


2  2  2 

a(cos     i)   -  2  cos  i   a{cos  i)   =  A(sin     i)  +  e 


;a8) 


Carrying  out  the  differentiation   indicated  by  the   symbol    a,   we  find  that 


Ai   =  -e  /(2   sin  i   cos   i) 


{A9) 


At  minimum  deviation,  a,  i,  and  e  are  related  by  the  equation 


+  a  =  2i, 


(AlO) 
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The  error  in  measuring  9  is  consequently  twice  Ai .  Likewise,  Fig.  A2  shows  that  a  levelling 
error  will  cause  a  similar  error  in  measuring  a. 

Suppose  that  the  prism  has  index  1.5  and  that  the  prism  angle  is  60°.   If  we  take  e  to 
be  0.1°,  then  we  find  from  Eq.  (A9)  that  Ai  is  about  0.1'  and  therefore  that  Aa  and  A9  are 
each  about  0.2'.  This  is  nearly  comparable  to  the  precision  of  the  instrument  used  in  the 
previous  examples.  Levelling  must  consequently  be  precise  to  better  than  0.1°. 

Levelling  error  is  a  systematic  error  and  is  always  negative  because  Ai  is  always 
negative.   It  results  in  an  underestimate  of  both  6  and  a,  and  therefore  in  an  underestimate 
of  n. 

The  index  of  a  liquid  may  be  measured  using  a  hollow  glass  prism.  The  measurement  will 
be  accurate  if  the  windows  of  the  prism  have  parallel  faces;  in  Appendix  B,  I  show  that  this 
requirement  is  not  stringent.  Therefore,  all  the  foregoing  remarks  refer  to  liquid-prism 
techniques  as  well  as  to  more-common  solid-prism  techniques. 
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Figure  A3.  Critical -angle  method  of 
determining  refractive  index  of  a 
liquid. 


It  is  also  possible  to  measure  the  index  of  a  liquid  using  a  prism  spectrometer  and  a 
critical-angle  method,  such  as  the  two  described  in  Ref.  54.  Both  methods  are  in  principle 
the  same;  the  geometry  is  illustrated  in  Fig.  A3.  The  liquid  to  be  measured  is  in  contact 
with  one  face  of  a  glass  prism  whose  index  is  known.  The  index  of  the  liquid  is  n  and  that 
of  the  glass  is  p.  u   must  be  greater  than  n  if  total  reflection  is  to  occur;  this  is  not 
generally  a  serious  constraint  where  fiber  optics  is  concerned,  because  the  liquids  used 
will  generally  have  indices  less  than  1.5. 

The  index  of  the  liquid  is  found  from  the  equation, 


» 


N  sin  a  -  sin  (ji  cos  a 


(All) 
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derived  in  Ref.    54.     6  is  the  critical    angle,   N^  =   ij2  _  g-jpZ  ^^   3p(j  other  parameters  are  as 
drawn   in  Fig.   A3. 

To  estimate  the  precision  of  the  critical-angle  methods,   we  differentiate  Eq.    (All), 
yielding  three  equations, 

An     -   (N  cos  a  +  sin   (j)  sin  a)    Aa,  (A12) 

a 

An  ^  =   -[(sin  a  sin  c])  cos   (}))/N  +  cos   a  cos   ^']   Acf),  (A13) 

if 

and 

An     =    (m/N)    sin  a  Au,  (A14) 

where  the  notation  is  as  before. 

If  we  use  the  same  example  as  above  (namely  that  a  =  60°  and  \i  -   1.5),  and  assume  that 
the  index  of  the  liquid  is  1.45,  we  find  that  6  is  about  75°  and  (f)  about  -23°.  Therefore, 

An  =  0.4  Aa, 
a 

An,  =  -0.2  A^,  (A15) 

and  An  =  0.9  Ap. 

If   Aa  =   0.15  mrad   (0.5'),    then  An     +  An     ~  0.0001.     We   learned   in  connection  with   Eq.    (A4) 

a  (p 

that  (J  could  be  determined  to  0.0002;  using  the  critical -angle  methods  adds  an  additional 
uncertainty  of  0.0001.  If  a  liquid  prism  is  available,  the  result  may  in  principle  be 
slightly  more  precise  than  a  determination  by  the  critical-angle  method.  However,  the 
liquid  prism  I  obtained  commercially  had  such  badly  warped  faces  that  I  decided  not  to  use 
it. 

Equations  (A2),  (A3),  (A12),  (A13)  and  (A14)  may  be  used  for  estimating  the  instru- 
mental limit  of  error  only.  This  represents  the  greatest  precision  that  can  be  expected 
from  a  particular  instrument.  A  given  measurement  will  have  to  be  carried  out  several  times 
before  the  statistical  limit  of  error  is  as  small  as  the  instrumental  limit. 

Finally,  the  index  of  typical  liquids  varies  considerably  with  temperature,  the  coef- 
ficient being  of  the  order  of  0.0004  K"-'-.  Therefore,  if  the  index  is  to  be  determined  to 
greater  precision  than  the  figures  assumed  in  this  appendix,  the  temperature  will  have  to  be 
controlled  to  within  a  few  tenths  of  a  kelvin. 
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APPENDIX   B.      USE  OF  THE   HOLLOW  PRISM 

A  hollow  prism  may  be  used  to  measure  the  index  of  refraction  of  a  liquid,   provided 
that  the  glass  windows  have  reasonably  parallel   faces.     To  determine  the  imprecision  that 
results  from  lack  of  parallelism,   we  assume  that  any  wedge  is  shared  equally  by  the  two 
windows  through  which  the  light  passes.     This  is,   of  course,   an  approximation,  but  it  allows 
us  to  draw  the  ray  through  the  liquid  parallel   to  the  base  of  the  prism  when  the  deviation 
is  a  minimum.     Only  half  the  prism  need  be  drawn,   as  in  Fig.   Bl.     In  that  figure,   the  angles 
are  defined  as  indicated.     The  angle  labeled  6  is  one-half  the  deviation  brought  about  by 
the  whole  prism. 


Figure   Bl .     Geometry  for  calculating  effect  of  wedge  in  faces  of  liquid  prism. 


If  the  window  surfaces  are  parallel,   then  the  angles  are  as  labeled  in  Fig.    Bl. 
Suppose  that  we  incline  the  inner  face  of  the  window  by  an  angle  Aa,  as  shown  by  the  dotted 
line.     Then  the  angles  change  from  lo  to   lo  +  AI3  and  so  on.     We  wish   to  calculate  the 
change  Ae  in  6  that  results  from  the  wedge  angle  Aa. 

To  begin,  we  apply  Snell's  law  to  the  refraction  at  the  liquid-glass  interface, 


n   sin   I^  =   M  sin   U 


(Bl) 


Differentiating  both  sides  of  Eq.    (Bl),  we  find  that 


Al'   =   (n  cos  lo/p  cos   U)   aK 


(B2) 


where  AI3  is  equal  to  Aa.  Similarly, 
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Al^  =  (u  COS  I2/COS  ip  Al^,  (B3) 


but  AI2  is  not  equal  to  Aa. 

We  learn  from  triangle  ABC  that 


and  therefore  that 


=13-1^.1^-12  (84) 


AG  =  AK  -  Al'  +  Al'  -  Al„  (B5) 

and 

Ae  =  [1  -  n  cos  Io/Im  cos  U)]^!^  +  (m  cos  I2/C0S  l'^   -  DaI^.      (36) 

Using  triangle  ACD,  we  find  that 

Al^  =  Al^  -  Aa,   (87) 

which  we  substitute  into  Eq.  (85).  Equation  (85)  becomes 

AG  =  [(n  cos  IJcos   U)  -  u](cos  Ip/cos  I')Aa  .  (88) 

Even  if  n  and  m  differ  by  0.1,  C0SI3'  is  very  close  to  COSI3,  so  Eq.  (88)  may  be  simplified 
to 

A9  =  (n  -  m)(cos  Ip/COS  I')Aa.  (B9) 

If  we  use  the   same  example  as  before,   and  let  n  -    u  be  0.05,   we  find  that 

2Ae  ~  0.15Aa  (810) 

where  2ao  is  the  total  error  of  measuring  the  deviation  and  where  Aa  is  the  wedge  angle  of 
each  window.  Common  optical  windows  will  have  wedge  of  20"  or  so;  thus,  2a9  is  about  3"  or 
0.05'.  Because  this  error  is  much  smaller  than  the  errors  noted  in  Appendix  A,  wedge  in  the 
windows  will  not  be  a  factor  unless  we  need  to  know  n  substantially  better  than  +0.0002. 
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APPENDIX  C.   RESOLUTION  LIMIT  AND  DEPTH  OF  FOCUS 

The  resolution  limit  of  a  microscope  objective  in  air  may  be  expressed  as  x/NA,  where 
NA  is  the  numerical  aperture  and  is  equal  to  sin  U;  U  is  the  angle  formed  by  the  marginal 
ray  and  the  axis,  as  shown  in  Fig.  CI  (see  Ref.  38,  Chapter  2). 

As  a  rule  (except  in  metallurgical  microscopes),  the  image  is  projected  into  a  medium 
whose  index  n  differs  from  that  of  air.  This  is  also  shown  in  Fig.  CI.  The  resolution 
limit  RL'  in  the  medium  is  equal  to  x'/sin  U' ,  where  x'  =  x/n  is  the  wavelength  in  the 
medium.  Using  Snell's  law  at  the  interface,  we  find  that  sin  U  =  n  sin  U'.  Therefore,  RL' 
-   RL.  The  change  of  U  is  precisely  compensated  by  the  wavelength  shift.  This  is  a  well- 
known  result,  namely  that  the  resolution  limit  is  x/NA,  independent  of  the  index  of  the 

surrounding  material. 

2        2 

m  The  depth  of  focus  DF  is  given  by  X/2(NA)  =  X/2sin  U.  (That  is,  the  image  is  in 

acceptable  focus  within  a  distance  ±DF  of  the  image  point.)  The  depth  of  focus  DF'  in  the 

2  2 

medium  is  therefore  x'/2  sin  U',  which  reduces  to  DF'  =  nDF  =  nx/2(NA)  .  The  change  of 

wavelength  does  not  compensate  the  change  of  angle,  and  depth  of  focus  is  larger  in  the 

medium  by  a  factor  of  n. 


\ 


Figure  CI.     Geometry  used  for  determining  the  resolution  limit  in  matching  fluid. 
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